Work fluctuation theorems for harmonic oscillators 
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The work fluctuations of an oscillator in contact with a thermostat and driven out of equilibrium 
by an external force are studied experimentally and theoretically within the context of Fluctuation 
Theorems (FTs). The oscillator dynamics is modeled by a second order Langevin equation. Both 
the transient and stationary state fluctuation theorems hold and the finite time corrections are very 
different from those of a first order Langevin equation. The periodic forcing of the oscillator is also 
studied; it presents new and unexpected short time convergences. Analytical expressions are given 
in all cases. 
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In this letter, we investigate, within the context of the 
Fluctuation Theorems (FTs), the work fluctuations of a 
harmonic oscillator in contact with a thermostat and 
driven out of equilibrium by an external force. First 
found in dynamical systems 111 and later extended to 
stochastic systems 0, S IE 0|i these conventional FTs 
give a relation between the probabilities to observe a pos- 
itive value of the (time averaged) "entropy production 
rate" and a negative one. This relation is of the form 
P{o~)/ P(—a) — exp[crr], where a and —a are equal but 
opposite values for the entropy production rate, P{o~) and 
P{— er) give their probabilities and r is the length of the 
interval over which a is measured. In these systems, the 
above mentioned FT is derived for a mathematical quan- 
tity <7, which has a form similar to that of the entropy 
production rate in Irreversible Thermodynamics Q. 

The proof of FTs is based on a certain number of hy- 
pothesis; experimenting on a real device is useful not only 
to check those hypothesis, but also to observe whether 
the predicted effects are observable or remain only a the- 
oretical tool. There are not many experimental tests of 
FTs. Some of them are performed in dynamical systems 
H in which the interpretation of the results is very dif- 
ficult. Other experiments are performed on stochastic 
systems, one on a Brownian particle in a moving optical 
trap and another on electrical circuits driven out of 
equilibrium by injecting in it a small current .10]. The 
last two systems are described by first order Langevin 
equations and the results agree with the predictions of 
ref.[EEl- A s f ar as we know no theoretical predictions 
are available for systems described by a second order 
Langevin equation. The test using an harmonic oscillator 
is particularly important because the harmonic oscillator 
is the basis of many physical processes. Indeed the gen- 
eral predictions of FTs are valid only for r — > oo and 
the corrections for finite r have been computed only for 
a first order Langevin dynamics. 

In the present letter, we address several important 
questions. We investigate first the Transient Fluctua- 
tion Theorem (TFT) of the total external work done 



on the system in the transient state, i.e., considering 
a time interval of duration t which starts immediately 
after the external force has been applied to the oscilla- 
tor. We then analyze the Stationary State Fluctuation 
Theorem (SSFT) which concerns fluctuations in the sta- 
tionary state, i.e., in intervals of duration t starting at a 
time long after the external force has been applied. We 
also study a new case of stationary behavior obtained 
when the system is driven periodically in time |ll| . In 
this case, which is actually a very important one, no the- 
oretical prediction is available. We show that the finite 
time corrections for SSFT are already very complex in 
both these rather simple situations. 

To test the FT we measure the out-of-equilibrium 
fluctuations of a harmonic oscillator whose damping is 
mainly produced by the viscosity of the surrounding fluid, 
which acts as a thermal bath of temperature T. We recall 
here only the main features of the experimental set-up, 
more details can be found in ref . [13, |l3| . The oscillator is 
a torsion pendulum composed by a brass wire and a glass 
mirror glued in the middle of this wire. It is enclosed in 
a cell filled by a water-glycerol solution at 60% concen- 
tration. The motion of this pendulum can be described 
by a second order Langevin equation: 

TJTI + v T77 + C0 = M + V^bTv, (1) 
d^ at 

where 9 is the angular displacement of the pendulum, I e ff 
is the total moment of inertia of the displaced masses, v 
is the oscillator viscous damping, C is the elastic tor- 
sional stiffness of the wire, M is the external torque, ks 
the Boltzmann constant and r\ the noise, delta-correlated 
in time. In our system the measured parameters are 
the stiffness C = 4.5 x 10 -4 Nmrad -1 , the resonant fre- 
quency f — y/C/I e g /(2tt) = 217Hz and the relaxation 
time r^T 1 = 2I e s/v = 9.5ms. The external torque M 
is applied by means of a tiny electric current J flowing 
in a coil glued behind the mirror. The coil is inside a 
static magnetic field, hence M oc J. The measurement 
of 9 is performed by a differential interferometer, which 
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FIG. 1: a) Typical driving torque applied to the oscillator; b) 
Response of the oscillator to the external torque (gray line). 
The dark line represents the mean response 9(t) to the applied 
torque M(t). 



uses two laser beams impinging on the pendulum mir- 
ror 0, 0| . The measurement noise is two orders of 
magnitude smaller than the thermal fluctuations of the 
pendulum. 8(t) is acquired with a resolution of 24 bits 
at a sampling rate of 8192Hz, which is about 40 times 
Jo- The calibration accuracy of the apparatus, tested at 
M = using the the Fluctuation Dissipation Theorem, 
is better than 3% (see 0). 

To study SSFT and TFT we apply to the oscillator a 
time dependent torque M(t) as depicted in Fig.^,, and 
we consider the work W T done by M(t) over a time r: 
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[M(t)-M(U)]-dt. (2) 



The TFT implies that U = whereas U > 3r Q for SSFT. 
As a second choice for M(t), the linear ramp with a rising 
time r r is replaced by a sinusoidal forcing; this leads to a 
new form of stationary state which has never been con- 
sidered in the context of FT. We examine first the linear 
forcing M(t) = M t/r r (Fig.QJi,)), with M a = 10.4 pN.m 
and T r — 0.1 s = 10.7 r a . The response of the oscilla- 
tor to this excitation is comparable to the thermal noise 
amplitude, as can be seen in Fig.JIJ)) where 6 it) is plot- 
ted during the same time interval of Fig.QJi)- Because 
of thermal noise the power W T injected into the system 
(cq|3J) is itself a strongly fluctuating quantity. 

We consider first the TFT. The probability density 
functions (PDF) P(W T ) of W T are plotted in Fig.Hi) for 
different values of r. We see that the PDF are Gaussian 
for all r and the mean value of W T is a few ksT. We 
also notice that the probability of having negative values 
of W is rather high for the small r. The function 
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is plotted in fig|2t>) • It is a linear function of W T for any 
r, that is S(W T ) = S(r) W T . Within experimental error, 
we measure the slope S(r) = 1. Thus for our harmonic 
oscillator the TFT is verified for any time r. This was 
expected 0, , and we give a derivation of this generic 
result for a second order Langevin dynamics at the end 
of the letter. 

We now consider the SSFT with U > 3r Q in eq[3 We 
find that the PDF of W T , plotted in Fig^), are Gaussian 
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FIG. 2: TFT. a) P(W T ) for TFT for various r/r a : 0.31 
(o), 1.015 (□), 2.09 (o) and 4.97 (x). Continuous lines are 
Gaussian fits, b) TFT; §(W r ) computed with the PDF of 
a). The straight continuous lines are fits with slope 1, i.e., 
£(t) = 1, Vr. 
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FIG. 3: SSFT with a ramp forcing, a) PDF of W T for var- 
ious r/r Q : 0.019 (o), 0.31 (□), 2.09 (o) and 4.97 (x). b) 
Corresponding functions S(W T ). c) The slope E(r) of S(W T ) 
is plotted versus r (□: experimental values; continuous line: 
theoretical prediction eq.QHJ with no adjustable parameter). 



with many negative values of W T for short r. The func- 
tion S(W T ), plotted in Figl^Js), is still a linear function of 
W T , but, in contrast to TFT, the slope £(t) depends on 
t. In FigEt) the measured values of S(r) are plotted as 
a function of r. The function S(r) — > 1 for t r a . Thus 
SSFT is verified only for large r. The finite time correc- 
tions of SSFT, which present oscillations whose frequency 
is close to fo, agree quite well with the theoretical pre- 
diction computed for a second order Langevin dynamics 
that we will discuss at the end of the paper (see ea llOjl . 
We stress that the finite time correction is in this case 
very different from that computed in ref. for the first 
order Langevin equation. 

The results of Figs l2l3l have been checked for several 
M /r r without noticing any difference. The errors bars 
in the figure are within the size of the symbols, and they 
come only from the calibration errors of the harmonic 
oscillator parameters, and statistics of realisations (typi- 
cally 5 x 10 5 cycles have been used). 

Finally we want to briefly describe the results of the 
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FIG. 4: Sinusoidal forcing, a) Sinusoidal driving torque ap- 
plied to the oscillator, b) Response of the oscillator to this 
periodic forcing (gray line) ; the dark line represents the mean 
response 9 it), c) PDF of the work W n integrated over n 
periods of forcing, with n = 7 (o), n = 15 (□), n — 25 
(o) and n = 50 (x). d) The function S(W n ) measured at 
uid/'iir = 64Hz is plotted as a function of W n for several n: 
(o)n = 7; (D)n = 15 (o)n = 25; (x)n = 50. Continuous lines 
are functions S(W„) computed from Gaussian fits of PDF (in 
Fig.b). e) The slopes E(n), plotted as a function of n for two 
different driving frequencies Ud = 64 Hz (□) and 256 Hz (o); 
continuous lines are theoretical predictions from ea llll with 
no adjustable parameter. 



periodic forcing. In this case M(t) — M sm.Wdt and the 
work expression (eq|2J) is replaced by 



W n = W T 
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M(t)-dt, (4) 



with r n = nflit/uid with n integer. This is a stationary 
state that has never been studied before in the context of 
FT. We find that for any driving frequency ui d the PDF 
of W n are Gaussian. The function S(W n ), measured at 
ujd/^ — 64Hz and plotted in Fig^Ji), is linear in W n 
and the corresponding slope S(rt) is a function of n. The 
measured values of £(n) are shown as function of n in 
figllJ}), where the results obtained at u>d/^ — 256Hz are 
plotted too. We see that the convergence rate is quite 
different in the two cases, which agree with our theo- 
retical predictions for a second order Langevin equation 
(see ea lll|) . Also in the case of the sinusoidal forcing the 
agreement between the computed and measured finite 
time corrections is very good. These results prove not 
only that FTs asymptotically hold for any kind of forc- 
ing, but also that finite time corrections strongly depend 
on the specific dynamics. In the case of the sinusoidal 
forcing, the convergence is very slow: in Fig.0^), we see 
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FIG. 5: a) PDF of the fluctuations 89 = 9 - 0(t) when the 
torque is applied (o), compared with a Gaussian fit of the 
PDF at equilibrium (continuous line), b) The measured spec- 
trum of 58 (o) is compared with the prediction of fluctuation 
dissipation theorem in equilibrium (continuous line). 



that it takes several dozens of excitation period (500 ms 
for LOdjl'K — 64Hz) to get S(n) = 1 within one percent. 
On the contrary for a ramp forcing, this was achieved 
after a few t q (20 ms), see Fig- Et) • 

Let us now compute the finite time corrections to TFT 
and SSFT (plotted in Figs|3];0^) for the harmonic oscil- 
lator, applying a method very similar to the one already 
used in the the context of the Jarzynski equality (see 
ref.0). We write 6{t) = 9(t) + 66(t) where 6{t) is the 
mean response of the system to the external torque and 
S9(t) are the thermal fluctuations. The mean response 
6{t) (dark line in Fig.^> resp. 2b) is computed by per- 
forming an ensemble average of 0(t) over 10 3 responses 
to the M(t) of Fig-Ht) resp. Fig.gji). It turns out that 
the measured 6 it) is equal to the solution of eq^ with 
X] = and with M equal to the applied time dependent 
torque. Once the mean behavior is known, it is useful 
to compare the statistical properties of 5 9 it) measured 
at M(t) 7^ with the equilibrium ones. In FigJSJi) we 
plot the gaussian fit of the PDF of 60 it) measured at 
equilibrium M(t) = (continuous line) and at M(t) ^ 
(o). The two curves are equal within experimental er- 
rors. Thus we conclude that the external driving does 
not perturb the equilibrium PDF, which is a Gaussian of 
variance UbT/C. In FigOJ>) we plot the power spectra 
of 56 in equilibrium (continuous line) and out equilib- 
rium (dotted line). The two spectra are equal and they 
coincide with the theoretical spectrum of an equilibrium 
second order Langevin dynamics (eq^with M — 0) com- 
puted using the oscillator parameters. Thus we clearly 
see that, within statistical accuracy, the fluctuations of 
89 measured at M(t) ^ are those of equilibrium. This 
important observation is the key point to estimate the 
finite time corrections of FTs. 

In order to compute S(r), we first decompose the total 
work of the external torque into the sum of a mean part 
and a fluctuating one, i.e. W T — W T + 8W Tl where 
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with a — 1 for the ramp forcing (eq. [3J or a = for si- 
nusoidal driving (eq. 0} . From those equations and the 
afore mentioned experimental observation on the fluctua- 
tions 89, we see that the fluctuations 5W T have a gaussian 
distribution, so the PDF of the work W T done by the ex- 
ternal torque is fully characterized by its mean W T and 
variance a 2 = ^(<5W r ) 2 ) where (.} stands for ensemble 
average. In such a case the FTs take a simple form: 



S(W T ) 



2W 



(7) 



where E(r) = (1 — e(r)) -1 following the notation of 
ref. 6] . From eq[7| it is clear that to estimate the finite 
time correction S(r) we need only to compute W T and 



a 2 . 



W T is simply obtained by inserting in eqJSJthe expres- 
sion of M(t) and the solution 9 of eq^with 77 = 0. The 
variance a 2 is computed using eq|S| For the linear ramp 
of figlD we obtain: 
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which can be computed using the correlation function 
R(r) = (S9(t)S9(0)). As already explained, the exper- 
imental data indicate that the statistical properties of 
59(t) on the ramp are the same as the properties of the 
equilibrium fluctuations, which are well described by a 
second order Langevin dynamics. Thus we can use for 
R(t) the known equilibrium correlation function of the 
thermal fluctuations which for a second order Langevin 
dynamics is [liij : 



R(r) 



k B T sin (ip\r\ + <-p) 
C shop 



exp(-a|r|), (9) 



where a = l/r a , a 2 + 'i/j 2 = u 2 = C/I c g and sinip = 
iP/loq. We checked our method on a first order Langevin 
equation, for which the exact results of refs.0, |g are 
available. We find that e(r) for the work computed with 
our technique in a first order Langevin dynamics is the 
same as in ref. HQ both for TFT and SSFT. Thus we 
can now safely apply our technique to a second order 
Langevin equation. We find that in the case of the TFT 
e = Vt, whereas in the case of the SSFT 
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The same calculations can be performed for any kind of 
M(t). For example with a sinusoidal forcing M(t) = 
M sm.(w d t), SSFT for the work W n defined in eq.^gives 



cos 27 U)\ 
2aT n 



(11) 



where 7 is the phase shift between 9(t) and M(t), i.e, 



tan(7) = — 2au>d/(u)Q — u>~j) and 



2mr/ujd with n 



integer. In Eq. ^2 0(e~ QT ") is a term that vanishes 
exponentially in crr n , the expression of which is compli- 
cated and will be reported in a longer article, together 
with many other interesting features. 

These analytical results agree remarkably well with the 
experimental results for TFT and SSFT for the work fluc- 
tuations in a harmonic oscillator (see Figs. 12: and^J?). 

In conclusion we have applied the FTs to the work 
fluctuations of an oscillator driven out of equilibrium by 
an external force. The TFT holds for any time whereas 
the SSFT presents a complex convergence to the asymp- 
totic behavior which strongly depends on the form of 
the driving. The exact formula of this convergence can 
be computed using several experimental evidences of the 
statistics of the fluctuation. These results are useful for 
many applications going from biological systems to nan- 
otechnology, where the harmonic oscillator is the simplest 
building block. 

This work has been partially supported by EEC con- 
tract DYGLAGEMM. 
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